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Abstract 

A variant of the one-dimensional t-J model with anisotropic spin interaction 
and broken parity is studied by the nested algebraic Bethe-ansatz method. 
The gapless charge excitations and the gapful spin excitations are obtained. 
It is shown that the singlet-superconducting correlation dominates in the low- 
density region by applying the finite-size scaling analysis in the conformal field 
theory. 



Typeset using REVT^ 



In these years, strongly-correlated electron systems have drawn much attention. It is 
partly caused by the discovery of the copper-oxide high-Tc superconductors. The one- 
dimensional models of strong correlation play an important role since some of them can 
be solved exactly. The non-perturbative results thus obtained contain some of the essential 
properties of the strongly correlated systems. Close scrutiny, however, is required when one 
tries to discuss higher dimensions based on the one- dimensional results, since some of them 
are specific to one dimension. 

Besides the Hubbard model the t-J model is regarded as one of the most 

basic models which contains the essence of strong correlation. The Hamiltonian is Titj = 

where n's are 



the number operators given by = + 71^ = cl^Ci^ + cj^Cii and the opm wpciatwio 
are S'f = | J2a,i3 cl^'^^,i3Cit^ with the usual Pauli matrices cr's. The Gutzwiller projector 



V = nj=i(l — '"'jT'^ii) restricts the Hilbert space by forbidding double occupancies hence 
represents strong correlation. In one dimension, the exact solution was obtained for the 
super symmetric case (2t = J) using the Bethe-ansatz method. The long-distance 

behavior of the correlation functions was also investigated by applying the finite-size scaling 
analysis in the conformal field theory to the excitation spectra obtained by the Bethe-ansatz 
method It was shown that the superconducting correlation does not exceed the others 
such as the spin density wave (SDW) and the charge density wave (CDW) correlations at any 
filling. However, the region where the superconducting correlation is dominant was found 
between the low-density supersymmetric region and the "phase-separated" region {J ^ t) 
by numerical diagonalization of finite clusters [|10| . Since the numerical results are not totally 
reliable due to the finite-size effect, it is highly desirable to have exact results. 

In this paper we study a one-parameter family of the one-dimensional correlated electron 
systems which includes the ordinary supersymmetric t-J model. The Hamiltonian is 



+7] (S^rii+i - ni^^i) + 2A n 



(1) 
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where — r^^ = 1 and we parameterize them as A = cosh 7 and 77 = sinh7 (7 G R>o) ||11|| . 
When 7 = 0, our Hamiltonian reduces to the ordinary supersymmetric t-J model. When 
the number of electrons coincides with the number of lattice sites, our model becomes s = ^ 
XXZ spin chain. The Hamiltonian (Q) commutes with the transfer matrix of the solvable 
two-dimensional classical lattice model associated with the supersymmetric quantum affine 
superalgebra f/g(s/(2|l)) (g = e"^), which is a special case of the model given by Perk and 
Schultz []T2| , p!3| . Thus one could say that the Hamiltonian is "g-supersymmetric" . The ground 
state, excitations and the correlation functions are obtained by the nested algebraic Bethe- 
ansatz method, the finite-size scaling analysis in the conformal field theory and the numerical 
diagonalizaton for small clusters. We find spin gap and observe that the superconducting 
correlation dominates in the low-density region in contrast to the ordinary supersymmetric 
t-J model. 

The third term of the Hamiltonian (|l]) breaks the parity invariance. The role of this term 
was investigated by diagonalizing the system without it numerically. It turned out that the 
spin gap survives in some cases. 



1. Bethe-Ansatz Equations (BAEs). The diagonalization of the Hamiltonian (P with 
periodic boundary condition reduces to solving the coupled algebraic equations (BAEs) 
derived by the nested algebraic Bethe-ansatz technique. They are 

/ sinfe + |7) \^ ^ 

JX. sin(Aa - Pj + f 7) 
j=i sin(A« - pj - f 7) 

where L is the number of lattice sites, N is the number of electrons, M is the number of 
down-electrons (magnons), p's are the quasi-momenta of electrons and A's are the magnon 
rapidities 0]. 

Hereafter we take the following ansatz for the ground state and the elementary exci- 
tations: A's are one-strings {Aq, G R|a = 1,---,M} and p's consist of the one-strings 
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TT sin(p,- - A^ -I- 17) 
i=i sin(Pi - A^ - |7) 

_ _ ^ sin(A^ - A^ + Z7) 
^J-^ sin(A„ - A/3 - 27) 



a 



1,2, 



1,2, 



(2) 
(3) 



{pj = Uj G R|j = 1, ■ ■ - 2M} and the two-strings {p^ = A„±27/2|a = 1, ■ ■ M}. Note 
that the real parts of the two-strings coincide with the magnon rapidities and, if M = N/2, 
there is no degrees of freedom for one strings of the quasi-momenta. This is essentially the 
same ansatz established in Refs. In our case, however, A's and m's are in the interval 

[— |, |] due to the periodicity of the BAEs (0), (^. By taking the logarithm of the BAEs, 
we have 

L i) = 2ml^ + E - A^, I) ^- = i, . . . , AT - 2M, (4) 

/3=1 

N-2M M 

L0(A„,7) = 27rzJ„+ ^ 0(A„ - m,-, ^) + ^ 0(Aa - A^, 7) a = l,---,M, (5) 

j=l ^ (3=1 

where 0(2;, a) = log (see ^^), and { Jj | j = 1, 2, ■ ■ ■ , — 2M} is a set of integers 

(or half-odd integers) if M is even (or odd) and the set { Jq, | a = 1, 2, ■ ■ ■ , M} is a set of 
integers (or half-odd integers) if A^ + M-|- 1 is even (or odd). We order the quantum numbers 
J's and J's according to Ij > Jj+i and Jq > Ja+i- 

In the thermodynamic limit L,N 00, the distributions of A's can be described by 
the continuous density given by L p(Aq,) = lim^^^Ar^^oo l/( Aq,+i — Aq, ). The energy and 
momentum up to 0(1) are 

E = I sinh7 E0'(P., |) = ^ sinh7 U j dk' p(A') 0'(A',7) + E |) I > (6) 

and 

Af „ N-2M 



P = iY.^{p,,l) = iL f d^ p{t^)cl>{K',^)+i Y^^i^vl)- (7) 
j=i ^ -I j=i ^ 

The integral intervals in (^ and (|^) will be discussed in the following sections. 



2. Ground State. We set A^ to be even for simplicity. Since the spin interaction is anti- 
ferromagnetic, the total 5*^ for the ground state is expected to be zero. This can be achieved 
by setting M = N/2, i.e. for the sector without one-strings pj = Uj. We also require that 
the momentum P to be zero. We propose the following ansatz for J's: the distributions 
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of J's for the ground state is restricted as Jmax > \Ja\ > Jmin, where Jmax = ^ ^ ^ o,i^d 



J If 



L-2M+1 



In the thermodynamic hmit, we assume that A's are distributed only in the regions 
[— 7r/2, ~Qg\ and [Qg, 7r/2] in accordance with the distribution of the quantum numbers J's. 
BAEs (H) and are reduced to 



2m p,(A) = -0'(A,7) + 



+ 

-Tv/2 



dA' p,(A') 0'(A-A',7), 



(8) 



where Qg is determined by /_^/2 + J^Qg Pgi-^) — ^'Hr = f • The ground-state energy E, 



N/2 



is given by 



Eg = i L sinh 7 



+ / 

-tt/2 J Qg . 



dA pg{A) 0'(A,7). 



(9) 



These equations can be solved numerically for arbitrary filling n. The results are given in 
Fig. 1. 

3. Charge Excitations. The charge excitations are those caused by the replacements 
of the A's while keeping M = N/2 namely remains to be zero. Thus the elementary 
excitations for the charge sector consists in making a jump (hole) at the point J^^ and putting 
a quantum number Jap^t a previously unoccupied region In the thermodynamic limit, 

BAEs @ and @ are reduced to 

Ott? 1 

2m pe(A) = -0'(A, 7) - — 5(A - Ah) + j<P\A - A„ 7) 



+ 



-tt/2 



+ 



■k/2- 



dA' p,(A') 0'(A-A',7), 



(10) 



retaining terms up to 0(L~^), where Ap and A^ denote the position of the hole and particle 
in the sea of two strings associated with the quantum number Jq,^ and Ja^ respectively. Qc 
is determined by /I^/2 + /qI^ Pc{A) = ^^'^^^"^ . For convenience, we decompose Pc(A) 
into the regular part and the singular part as Pc(A) = Pco(A) — -^pci(A) — -^^{A — A^), where 
Pco(A) satisfies 

-Qc r7T/2-\ 

+ / 

-n/2 J < 
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2m pco(A) = -0'(A,7) + 



dA' p,o(A') 0'(A-A',7). 



(11) 



For Pel (A), we have 
2m Pel (A) = 0'(A - A;,, 7) - 0'(A - Ap, 7) + 



-tt/2 



+ 



n/2- 



rfA'pei(A')0'(A-A',7) (12) 



The excitation energy AE from the ground state and the momentum P are given by 

rn/2' 



AE = i sinh7 |^0'(Ap,7) - 0'(A;„7) 



P = i U(Ap,7)-0(A,,7) 



-tt/2 



-tt/2 



+ 



tt/2 



rfApei(A) 0'(A,7) 



rfApei(A) 0(A,7) . 



(13) 



(14) 



Solving ([I2D numerically, the dispersion relation for the elementary charge excitations was 
obtained. The result for 7 = 2 and n = 0.45 is shown in Fig. 2. The results for other 
parameters do not change in an essential manner, namely the charge excitation is always 
gapless. 



4. Spin Excitations. The spin excitations can be considered as excitations coming from 
destroying the two-strings p^'s and creating one-strings m's. To study the elementary ones, 
let us consider the case of M = N/2 — 1 magnons. We assume that, in the sea of the quantum 
numbers J's there are no jumps 0J^. Then, in the thermodynamic limit, BAE (§) becomes 



2m p,(A) = -0'(A, 7) + ^0'(A - «i, |) + ^0'(A - ^2, |) 



+ 



-tt/2 



+ 



tt/2' 



dA' p,(A') 0'(A-A',7), 



(15) 



where ui and U2 are one-string quasi-momenta, and Qs is determined by 
1-^/2 + Iq^ dAps{A) = It is convenient to decompose p<i(A) into contributions 

of order 0(1) and 0{L~^) as Ps(A) = pso(A) — -^Psi(A), where Pso(A) satisfies the same 
equation as ([TT|) obtained by replacing all the suffices c to s. Then the integral equation for 
Psi(A) is obtained as 



2m p,i(A) = -0'(A - m, |) - 0'(A - U2, |) + 



-Qs rTT/2' 
+ 

-tt/2 J Qs . 



rfA'p,i(A')0'(A-A',7), 



(16) 



The excitation energy AE" from the ground state and the momentum P are 

l"!r/2' 



AE = i sinh7 



-n/2 



+ 



7 7 



+ 



Tr/2 



rfAp,i(A)0'(A,7) (17) 



rfAp,i(A)0(A,7) . (18) 



Solving (|T6|) numerically, the dispersion relation for the elementary spin excitations was 
obtained. The result for 7 = 2 and n = 0.45 is shown in Fig. 3. The results for other 
parameters do not change in an essential manner, namely the spin excitation is always 
gapful. The spin gap as a function of 7 is also shown in Fig. 4. and one can see that the 
gap increases as holes are doped. 

To study the effect of the parity violating term, we calculated the spin-spin correlations 
{S-Sj) for the parity-unbroken Hamiltonian Htj — J^iVi^i^i+i ~ ^i^i+i) 1131 • "^^^ results 
are shown in Fig. 5 and they indicate that the correlation decays exponentially. Hence there 
is still excitation gap. 



5. Correlation Functions. Consider a field-theoretic description of the low-lying exci- 
tations. Since the dispersion for the low energy charge sector is approximately linear for 
< n < 1, and the gapful spin sector is irrelevant for the low-energy behavior, we can expect 
the system can be described by the conformal field theory |[T8| . 

Let us consider the excitations described by the density p(A) satisfying 



2TTip{A) = -0'(A,7) + 



f"^- + r^']dA'p(A')0'(A-A',7), (19) 



and apply the general method of Kawakami-Yang p| for the finite-size scaling method |119| , |20 
Using the Fourier-transform technique, we rewrite (0) as 

p(A) = 2i?,(2A) + f^^dA' 2i?g(2(A - A'))p(A'), (20) 

•/ Q — 

where we have introduced the deformed Shiba-function defined by Rq{v) = 

^ Emez i+J^ i li • The energy is given by 



E/L = 2 cosh 7 — 27r sinh7 



2i?,(0) 



dK 2Rq{-2A) p(A) 



(21) 



Thanks to (pOf ) and (0), we can immediately apply the general argument and the results 
are: i) the charge sector can be described by the c = 1 bosonic conformal field theory 
i.e. it belongs to the universality class called the Tomonaga-Luttinger liquid, ii) the com- 
pactification radius is given by r = ^{Q), where the dressed charge ^{K) satisfies 

^(A) = l+S%dk' 2i?g(2(A-A')) ^(r^), and Q is determined by [/:J/2 + /q ^] p(A) = 

As usual, we parameterize r by Kp = r^/2. The equation for ^ was solved numerically 
and Kp as functions of n are shown in Fig. 6. The relations between Kp and the critical 
exponents are shown in Table I As long as 7 7^ 0, the singlet-superconducting 

correlation is dominant when Kp > 1 in the low density region (high doping). However 
for the usual supersymmetric case (7 = 0), the superconducting correlation can not be 
dominant in any filling as seen in Fig. 6. This indicates that deformed t-J models 
including ours may be more appropriate to study superconducting mechanisms than the 
ordinary t-J model. 

The authors benefited from the help in the numerical calculations by Y. Hatsugai and 
Keiko Kimura. J. S. thanks Kazuhiro Kimura for discussions on the deformed t-J model 
Hamiltonian. They are also grateful to M. Ogata and S.-K. Yang for many useful comments 
and suggestions on the algebraic Bethe ansatz and the finite-size scaling analysis. 
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FIGURES 



Fig. 1. Ground state energy per site as a function of electron density n. The solid lines 
are obtained by solving We also performed direct numerical diagonalizations of the 
Hamiltonian for L = 14 with 7 = 2.0, 1.0 and 0.5. The results are plotted by □ , + and o 
respectively. 

Fig. 2. Dispersion of the elementary excitations in the charge sector for 7 = 2 and 
n = 0.45. Sufficiently many points in the continuous spectrum are shown. The momentum 
P is periodic with period 27r. The gapless points are at P = 0, 2kF and 27r — 2kF- 

Fig. 3. Dispersion of the elementary excitations in the spin sector for 7 = 2 and 
n = 0.45. Sufficiently many points in the continuous spectrum are shown. The momentum 
P is periodic with period 2tt. There is no gapless point. Note that the momentum P for the 
elementary spin excitation is restrected in [—Pm,Pm], where Pm depends on 7 and n. 

Fig. 4. Spin gap as a function of 7. For n = 1, the spin gap of our model reduces to 
that of the XXZ spin chain (see Eq.(g)) [g. 

Fig. 5. Spin-spin correlations on a logarithmic scale for the parity-unbroken Hamiltonian 
'Htj — J2iV{Si''^i+i~^iSi+i)- The results are obtained by numerical diagonalization of L = 12 
systems for (a) = 10, 7 = 1.5, (b) = 8, 7 = 1.5 and (c) = 2, 7 = 2.0. 

Fig. 6 Kp{n)^s are shown for 7 = 0.0, 0.5, 1.0 and 2.0. The broken line (for 7 = i.e. the 
ordinary supersymmetric case) denotes the data from Ref. 8. It can be shown analytically 
that Kp{0) = 2 and Kp{l) = 1/2 for any 7 > 0. 
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TABLES 

TABLE I. Relation between Kp and the critical expo- 
nents of the correlation functions. 



Correlations 



exponents 



2kF SDW (spin density wave) 
2kF CDW (charge density wave) 
SS (singlet superconductivity) 
TS (triplet superconductivity) 
4kF CDW (charge density wave) 



exponential decay 
exponential decay 
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